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^ . ABSTRACT 

C/3 _ 

I Context. There is growing evidence of relativistic jets in radio-loud narrow-line Seyfert 1 (RL-NLS 1) galaxies. 

Aims. We constrain the observational properties of the radio emission in the first RL-NLS 1 galaxy ever detected in gamma-rays, PMN 
J0948+0022, i.e., its flux density and structure in total intensity and in polarization, its compactness, and variability. 
Methods. We performed three real-time e-VLBI observations of PMN J0948-F0022 at 22 GHz, using a global array including tele- 
scopes in Europe, East Asia, and Australia. These are the first e-VLBI science observations ever carried out with a global array, 
reaching a maximum baseline length of 12458 km. The observations were part of a large multiwavelength campaign in 2009. 
Results. The source is detected at all three epochs. The structure is dominated by a bright component, more compact than 55 pas, 
with a fainter component at a position angle 9 ~ 35°. Relativistic beaming is required by the observed brightness temperature of 
3.4 X 10" K. Polarization is detected at a level of about 1%. 
CNj , Conclusions. The parameters derived by the VLBI observations, in addition to the broad-band properties, confirm that PMN 

c 2 ^ . J0948+0022 is similar to flat spectrum radio quasars. Global e-VLBI is a reliable and promising technique for future studies. 



X 



Key words. Galaxies; active, Seyfert, jets - Instrumentation: high angular resolution, interferometers 

1. Introduction relativ istic jets (IZhou et al.l2003l;lkomossa et al.l2006tlDoi et al.l 

l2006h . 

The detection by Fermi/LAT of gamma-rays from a handful 

Radio-loud (RL) narrow-line Seyfert 1 (NLSl) active galactic of RL-NLS 1 has provided definitive confirma tion of the pre s- 

nuclei (AGNs) have received considerable attention in recent ence of relativistic je ts in these sources (Abd o etal.ir2009blldL 

years, since they pose a challenge to current unified schemes iFoschini et al.l 1201 Oah . A multiwavelength (MWL) campaign 

(iUrrv&Padovanilll995l) . NLSl AGNs are characterized by an was organized in March-July 2009 to investigate the properties 

optical spectrum with narrow permitted lines FWHM(Hj6) < of PMN J0948-I-0022 with simultaneous observations covering 

2000 km/s, a ratio of the equivalent widths of [OIII]/15007 to Hfi gamma-rays. X-rays, UV, optical, and radio. This campaign has 

smaller than 3, and a bump caused by Fell (see, e.g.. |Poggd2000l delivered some intere sting results on th e emission processes in 

for a rev iew). They also exhibi t prominent soft X-ray excesses this z - 0.585 source (lAbdo et al.l2009d ). The direct observation 

jBoUer et al. 1996; Grupel l2004l) . These properties are indicative of simultaneous variability across the electromagnetic spectrum 

of very high (near Eddington) accretion rates and relatively low has confirmed the identification of the gamma-ray source with 

black hole masses (10^ - 10*^ Mq). Only a small percentage (~ the lower energy counterpart. This coordinated variability is also 

7%) of NLS 1 are radio-loud. In these cases, their flat radio spec- characteristic of the relativistic jets present in the central regions 

tra and VLBI variability suggest that several of them could host of blazars. Overall, the multiwavelength observational proper- 
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Table 1. Log of observations. The two-letter station abbreviations are: AT - Australia Telescope Compact Array (ATCA, 6x22m), 
Cm - Cambridge (32m, limited to 128 Mbps in all experiments), Da - Damhall (25m, limited to 128 Mbps in all experiments), 
Ef - EfFelsberg (100m), Ho - Hobart (26m, limited to 128 Mbps and a single polarization in all experiments), Jbl - Jodrell Bank 
(Lovell Telescope, 76m), Jb2 - Jodrell Bank (Markll telescope, 32x25m), Ks - Kashima (34m, limited to 128 Mbps and a single 
polarization in all experiments), Mc - Medicina (32m), Mh-Metsahovi (14m), Mp -Mopra (22m), On - Onsala (20m), Pa - Parkes 
(64m), Sh - Shanghai (25m), Tr - Torun (32m), Ys - Yebes (40m), Wb - Westerbork Synthesis Radio Telescope (WSRT, 12x25m) 



Epoch # 


Day of 


Freq. 


Participating telescopes 


Longest baseline 


beam FWHM 


S 






2009 


(GHz) 




(km) 


(mas X mas, °) 


(mJy) 







Apr 21 


1.6 


Ef, Cm(+Da), Jbl, Mc, On, Tr, Wb 


Jb-Tr, 1 388 


35.4x22.9,11.7 


173 ± 


20 


1 


May 23 


22.2 


Sh, Mp, Ho, Ks, Ef, Mc, On, Cm, Jb2 


Ef-Ho, 12 359 


0.38 x0.15,39.1 


665 ± 


130 


2 


Jun 10 


22.2 


Sh, Mp, Ho, Ks, Ef, Mc, On, Mh, Ys, Cm, Jb2 


Ys-Mp, 12458 


0.43 X 0.14, 43.3 


345 ± 


70 


3 


Jul 04 


22.2 


Sh, Ks, Mp, Ho, Pa, AT, Ef, Mc, On, Mh, Cm, Jb2 


Ks-Mc, 8 811 


0.42 x0.35,39.5 


445 ± 


90 
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Fig. 1. Visibility amplitude versus (M,v)-radius for the observations on 2010 May 23. Stokes 'LL' is shown, since Hobart only 
observed at LCP The red points show the amplitude visibility of a Gaussian component with flux S = 665 mJy and FWHM 
0.205 x 0.055 mas in PA 53.1°. 



ties of J0948+0022 seem therefore to indicate that the physics 
of RL-NLS 1 and blazars are similar, and in particular a relation- 
ship between RL-NLS 1 and the powerful, high accretion, flat- 
spectrum radio quasars. 

In this letter, we present additional details of the very long 
baseline interferometry (VLBI) observations carried out during 
this campaign. Section 2 describes the procedures of the real- 
time e-VLBI observations (ISzom oru 2008), which were the first 
application of this technique using a global array of antennas 
in an astronomical program. Section 3 presents the results, and 
Sect. 4 discusses conclusions and prospects. 

2. Observations and data reduction 

The observations took place on 2009 April 21, May 23, June 10, 
and July 4 using a subset of the European VLBI Network (EVN) 
and the Long Baseline Array (LBA) in Australia, and the 34 m 
Kashima telescope in Japan. For a complete list of participat- 
ing telescopes and other observational parameters, we refer to 
Table [1] The initial observations on 21 April were performed at 
1 .6 GHz for two hours with the EVN only at 5 12 Mbps data rate. 
These observations served as a pilot project designed to assess 
the feasibility of the experiment; we refer to them as epoch 0. 

The other three experiments (epochs 1,2,3) were carried out 
with a global array at 22 GHz for maximum resolution, although 
Australia-Europe baselines could not be formed at the third 
epoch because of the availability constraints of the Australian 
telescopes. Each observation lasted for about 1 1 hours and was 



preceded by a 2-3 hour preparation and clock-search run for 
the Australian- Asian array, and later included an hour of clock- 
search and re-referencing for the whole array including the EVN. 
The fringe-finders used were 0420-014, OJ287, and 1055+018. 
Additional nearby calibrators 0736H-017 and 0805-077 were 
regularly observed for continuous fringe monitoring; these later 
also served as amplitude calibration check sources. The total 
flux densities of these compact sources were obtained from 
single-dish EfFelsberg (2nd epoch) and synthesis-array ATCA 
(3rd epoch) measurements. The total bit rate per telescope was 
512 Mbps, divided into four 16 MHz sub-bands in both polar- 
izations. The data were stre amed in real-time fro m the stations 
to the EVN data processor (ISchilizzi et al.ll2001l) at JIVE in all 
experiments. Plots of the real-time frin ges as well a s the final 
coverage of the {u, v)-plane are shown in lGiroletti et a l. (2010c!). 
The longest baseline in our array was the Yebes-Mopra on June 
10, reaching 12458 km. 

The data were processed using the NRAO Astronomical 
Image Processing System (AIPS). Amplitude calibration was 
done with the measured gain curves and system temperatures. 
For the LBA stations, nominal SEFDs were used. The Tjys 
data were corrected for atmospheric opacity assuming standard 
weather conditions for a subset of stations that had opacity-free 
calibration data. This in general worked well except for the low- 
est elevations below 10°. All the sources were fringe-fit sepa- 
rately, no phase-referencing being applied. All of the stations 
produced fringes, but Kashima had setup problems in the first 
two experiments and Hobart had a setup problem in the first 
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Fig. 2. Total intensity contours for J0948+0022 at epoch 

2, with levels traced at (-1,1,2,4,...) x l.SmJybeam"'. 
The central pixel corresponds to RA= 09H8"57.320', 
Dec=+0° 22' 25.560". 

part of the first experiment. After fringe-fitting, the station am- 
plitudes were further adjusted using the calibrator sources. Phase 
and amplitude self-calibration were performed in Difmap to pro- 
duce the final images. 

The phase response of the right- and left-handed polarized 
signals with parallactic angle was corrected before fringe -fitting 
for all stations, including the new 40m Yebes radio telescope. 
For this Nasmyth offset mount antenna, we used the routines 
of the 31DEC1Q version of AIPS, following the procedures de- 
scribed in Dodson (2009); no RR/LL phase difference was seen 
on Yebes baselines after this calibration. The observations had 
not originally been designed to analyze full cross-polarizations; 
however, we tried to solve for the D-terms using the bright cali- 
brator 1055-1-018, limiting ourselves to the European stations at 
the second epoch (when there were the largest number of sta- 
tions). Reasonable solutions were found, which allowed us to 
study the source polarized flux. We did not however correct for 
the polarization angle because of the unknown R-L phase differ- 
ence at the reference antenna. 

3. Results 

The source is clearly detected at all epochs and frequencies. 
Inspection of both the visibility plots and the image plane shows 
that the structure is dominated by a bright compact component. 
The flux density is difficult to constrain to better than 20% be- 
cause of a somewhat incomplete information about the system 
temperature and gain curve from some stations. It is clear how- 
ever that there is a flux density of several hundreds milliJanskys 
revealed at all epochs, with possibly some variation. The over- 
all non-simultaneous spectral index between 1 .6 and 22 GHz is 
clearly inverted, with a = 0.3 - 0.5 (5(v) oc v"). Flux densities 
from visibility model fitting at all epochs are given in Col. 7 of 
Table □ 

In Fig. [U we plot the visibility amplitude as a function of 
(m, v')-distance at the first epoch, which is the one with the longest 
mutual visibility on the longest baselines. Despite the uncer- 
tainty in the absolute flux level, the relative amplitude scales 
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Fig. 3. Polarization image from epoch 2, European baselines 
only. Green contours show total intensity at (-1, 1,2,4,...) x 
1.2mJybeam"', with a peak of 250 mJy beam"'; white vectors 
show total polarization (1 mas = 1.2 mJy beam"', PA not cor- 
rected); color shades show fractional polarization (between 
and 1.5%). 



between telescopes in Europe and Australia are reliably cali- 
brated, as shown by the consistent values across the short (m, v)- 
distances. We are therefore confident that the ratio of the visibil- 
ity amplitudes on the longest to short baselines is reliably mea- 
sured, which suggests that that some substructure is present. The 
fitting of an elliptical Gaussian model to the visibilities yields 
a component size of 0.205 x 0.055 mas at PA53.1°. The minor 
axis corresponds to about one third of the synthesized beam; this 
means that the source was indeed resolved, thanks to the length 
of baselines in this intercontinental observation (up to 12 359 km 
between Effelsberg in Europe and Hobart in Australia). 

Moreover, the baselines to Shanghai provide visibilities at 
an intermediate range (~ 500 - 600 MA) in NS and EW direc- 
tion. A smaller visibility amplitude relative to that for the longest 
baselines indicates that there is structure in the direction trans- 
verse to the longest baselines (PA ~ -45°). In the second epoch, 
the Fourier transform of the visibility plane, cleaning, and self- 
calibration are able to more clearly constrain the total intensity 
structure. A jet component is found at r = 0.86 mas from the 
core in PA 6* = 35.4° (see Fig.|2]i. This is also consistent with the 
lower frequen cy structure, which shows a resolved jet feature in 
PA ~ 3 0-40° (iDoi et alJ2006l:lAbdo et al.l2009cl:lFoschini et al.1 
l20T0bh . 

The brightness temperature in the source rest-frame based on 
the model-fit is ~ 3.4 x 10" K for the first epoch, whi ch is com- 
parable to the equipartition brightness temperature (iReadheadI 
1994). Given the length of the minor axis derived by the model 
fit, it is conceivable that the size of the core along the major 
axis is not well-constrained only because of the relatively insufli- 
cient baseline lengths. Therefore, the major axis of the elliptical- 
Gaussian model component represents only an upper limit, and 
the resulting brightness temperature is only a lower limit. This 
confirms that the source must be relativistically beamed, since 
the obse rved T\, exceeds the equipartition brightness tempera- 
ture (see lDoi et al.]l2006HZhou et al.ll2003h . 
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Polarization could only be studied for the second epoch and 
using intra-European baselines. Figure|3]provides a polarization 
map for this epoch. The mean fractional polarization is 0.9% 
(peak of 1.3%), which is barely in excess of the one reported 
from MOJAVE at 15 GHz at the same epoch. The absolute polar- 
ization angle cannot be determined for our data. To our knowl- 
edge, this is the first report of a VLBI polarization experiment 
using the new 40m Yebes radio telescope. 

4. Discussion and conclusions 

4.1. Scientific impiications of the resuits 

The census of extragalactic gamma-ray sources throughout the 
EGRET era has been dominated by blazars (flat spectrum radio 
quasars and BL Lac type objects). The discovery of gamma-ray 
emission from radio-loud narrow-line Seyfertl nuclei by Fermi 
is therefo re of great importan ce to the study of relativistic jets 
in AGNs dAbdo et al.ll2009dh . First, it has confirmed that these 
relativist ic structures exist in radio loud NLSl, as proposed by 
IZhou et al. (2003) and Doi et al. (2006). Moreover, exploiting 
the sensitivity and the surveying capability of Fermi, it has trig- 
gered the organization of the large MWL ca mpaign including th e 
global e-VLBI observations presented here (lAbdo et al.ll2009d) . 

The MWL campaign has revealed variability at all energies, 
and permitted the estimate of the physical parameters of the jet, 
such as the dissipation radius (r - 6.7 x 10'^ cm ), the mag- 
netic fi eld (B - 4.1 G), and the jet total power (lAbdo et all 
I2009ch . A bulk Lorentz factor F > 1 was implied by the large 
Compton dominance; in particular, F = 10 was assumed as a 
typical value. Our lower limit to the brightness temperatures 
estimated by elliptical Gaussian model fitting to the global e- 
VLBI visibihties is a few xlO" K and provides an independent 
confirma tion that this jet has a Doppler fac tor greater than one 
(see also Doi et al. 2006; Zhou et ani2003h . as well as in other 
RL-NLSl (G u & Ch en 2010). Because of the uncertainty in the 
overall amplitude scaling, a Doppler factor estimate based on the 
variability brightness temperature is not discussed in this study. 
However, a trend in the total flux density from our observation is 
present, with a peak measured on 2009 May 23 and a decrease 
in the following epochs, which are also characteristic features 
of blazar-like emission. This result also agrees with other inde- 
pendent variability measurements requiring S > 1 reported at 
various frequencies in Abdo et al. (2009c). 

Significant fractional polarization is also a characteristic of 
Fermi detected lets tHovattaet al.. 2010). In J0948+0022 itself, 
iFoschini et al.l ( 1201 Obi) reported an increase in the fractional po- 
larization up to 3% and a swing of the EVPA of about 90° in 
possible connection with the 2010 gamma-ray outburst. This 
is similar to what has been repo rted for other blazars, such as 
PKS 1502+106 ( Abdo et al.ll20r0 ). and might be related to the 
physical mechanisms responsible for enhanced gamma-ray ac- 
tivity. In this context, our detected value of ~ 1% fractional 
polarization seems more characteristic of an intermediate state 
of activity, consistent with the transition to a low state at high 
energy ob served towards the end of the 2009 MWL campaign 
(lAbdo et al. 2009c). 

4.2. Prospects for tiie use of global e-VLBI 

In this paper, we have presented radio interferometric observa- 
tions in real-time e-VLBI mode (Szomoru 2008), using, for the 
first time, a global array spanning three continents. The feasibil- 
ity of transferting data across continents and correlating them in 



real time has been demonstrated in a number of tests and in par- 
ticular in the LBA observations of SN 1987A, which were cor- 
related at JIVE (Tinga y et a l. 2009; Ng et"aL| l201 lb . However, 
these are the first astronomical e-VLBI science observations on 
a global scale. 

From the very beginning of VLBI, global baselines have 
been routinely achieved between e.g., Europe and the US, and 
some examples exist of truly global coll aborations between ar - 
rays on three different continents (see e.g. lFomalont et al.l2001h . 
However, it has always been a challenge to organise these ex- 
periments, especially in the case of time-critical observations 
such as our coordinated MWL campaign. The advantages of per- 
forming such observations in real time e-VLBI (when transport- 
ing the data through high-speed networks is possible from all 
participating telescopes) include the ability to monitor the net- 
work performa nce in real-time and th e possibility of delivering 
prompt results (iGiroletti et al.l l2010a b). But most importantly, 
with the standardiza^tion o f the data formats and transport proto- 
cols (IWhitnev et alJl2009l) . it will become much easier to create 
ad-hoc global networks at relatively short notice, providing ex- 
cellent Mv-coverage and unique long baselines, as well as greatly 
increased sensitivity due to the increase in sustainable data rates 
and the number of telescopes. 

Our pioneering global e-VLBI efforts, in spite of some initial 
technical difficulties, have convincingly demonstrated the oper- 
ational feasibility of this type of observation. Using global base- 
lines and a high observing frequency, we have probed the inner- 
most region at or above the self-absorption frequency, thus ob- 
taining the most valuable information as far as the comparison 
to high energy activity is concerned. The overall success of the 
MWL campaign and in particular of the global e-VLBI observa- 
tions presented here is highly encouraging for the continuation 
of the synergy between high energy astrophysics and VLBI. 
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